Atom detection in a two-mode optical cavity with intermediate coupling: 

Autocorrelation studies 
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We use an optical cavity in the regime of intermediate coupling between atom and cavity mode 
to detect single moving atoms. Degenerate polarization modes allow excitation of the atoms in one 
mode and collection of spontaneous emission in the other, while keeping separate the two sources of 
light; we obtain a higher confidence and efficiency of detection by adding cavity-enhanced Faraday 
rotation. Both methods greatly benefit from coincidence detection of photons, attaining fidelities 
in excess of 99 % in less than f fis. Detailed studies of the second-order intensity autocorrelation 
function of light from the signal mode reveal evidence of antibunched photon emissions and the 
dynamics of single-atom transits. 
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I. INTRODUCTION 

Atoms coupled to a finite number of modes of a cav- 
ity form a cavity quantum electro dynamic (QED) system 
These objects have intrinsic interest as physical sys- 
tems, and they are important for applications in quantum 
information science 0, H, 0, H[ . Beyond quantum infor- 
mation, single atom detection with cavities is the tool of 
choice for the study of atom dynamics and coherent atom 

optics @, 0,1, Hi Eg. 

Cavities have also enabled the study of quantum optics 
effects difficult to observe in free space. For example, Ref. 
[H| shows a transition from antibunching to bunching in 
light collected from an atomic beam into a cavity mode. 
The effect is visible because the cavity supports only one 
spatial mode out of the infinity in free space, and this 
ensures that the light detected from different atoms is 
spatially coherent. 

Any measurement protocol involving individual atoms 
in a cavity requires information that the atom is cou- 
pled to the cavity mode. The information manifests it- 
self in two ways: as an increase in the amount of light 
in the mode owing to collection of atomic spontaneous 
emission, or as a decrease in the amount of light that 
the cavity transmits, owing to atomic absorption or dis- 
persion. Cavity QED enhances the rate of spontaneous 
emission collection, and excitation perpendicular to the 
cavity mode enables high rejection of background light 
EH EH 13- Detection of a decrease in driven cavity 
transmission requires the averaging of the intrinsic shot 
noise until the intensity level change is significant, and 
needs strong atom-cavity coupling for rapid detection 
flU [l6l U3, EH ■ Experiments have demonstrated single 
atom detection using both methods for moving atoms 
d,El and trapped atoms EI El- 

We have developed a technique for atom detection that 
combines elements of driven cavity QED detection and 
spontaneous emission collection. We observe atoms in 
the mode volume of an optical cavity through the collec- 
tion of spontaneous emission and Faraday-rotated for- 



ward scattering into an orthogonally polarized cavity 
mode that is not driven Q • The technique allows rapid 
identification of highly-coupled atoms, as well as mea- 
surements of photon correlations and single-atom dynam- 
ics that are difficult to perform in free space. This paper 
presents our approach, in particular detailed autocorre- 
lation studies of the transmitted light that is basically 
the result of resonance fluorescence. We identify contri- 
butions from different physical processes to the autocor- 
relation signal in the short time (atomic decay) and long 
time (transit of atoms through the mode) [2l|, [22, [23| . 

The paper is divided as follows: Section [TTJ introduces 
the system, with a simplified theoretical model in section 
IIIII and a description of the apparatus in section [TV] Ex- 
amples of our raw signals follow in section [V] Section 
IVII explores the parameter space to establish the con- 
figuration for optimal fidelity (least chance of error in 
detection) . Section I VIII presents a study of the autocor- 
relation function of the perpendicular (undriven) mode. 
We conclude in section [Villi 



II. DETECTION SCHEME 

The basic components of our system are typical of cav- 
ity QED setups (see Fig.[T|). We couple a linearly polar- 
ized laser (the drive) to a TEMoo mode of a Fabry-Perot 
optical cavity resonant with the Z?2 line in 85 Rb. In- 
tersecting the cavity axis at near-normal incidence is a 
slow-moving beam of neutral 85 Rb atoms, each of which 
interacts with the cavity mode for some time before leav- 
ing. Because of the large difference in time scales (5 
/is transit time versus 26 ns excited state lifetime), the 
atoms become excited and scatter the drive light many 
times while in the mode volume. 

The crux of our detection scheme is the separation of 
the cavity output field into two distinct components rel- 
ative to the polarization of the exciting laser: one paral- 
lel (driven mode) which is little affected by the transit- 
ing atoms, and the other perpendicular (undriven mode) 
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FIG. 1: (Color online). Schematic of the apparatus with 
the basic optical elements necessary for coincidence measure- 
ments. Vertically polarized laser light (parallel) drives the 
high finesse cavity traversed by a continuous beam of cold 
85 Rb atoms from a magneto-optical trap. The output light 
passes through a polarizing beam splitter (PBS) that sepa- 
rates the horizontal (perpendicular) polarization, sending it to 
a nonpolarizing beam splitter (BS) for coincidence measure- 
ments using two avalanche photodiodes (APD). The direction 
of the magnetic field (B) inside the cavity can be parallel to V 
(7r-polarized drive) or parallel to the direction of propagation 
(a-polarized drive). 




m = n-1 n n+1 m = n-1 n n+1 
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FIG. 2: (Color online), a) A linearly polarized driving field 
excites the atom, but a transition to a ground state level with 
different m results in emission in the orthogonal polarization, 
b) Simplified diagram of 85 Rb levels relevant to Faraday ro- 
tation in an applied magnetic field. 



which in the ideal case is populated only by scattered 
light from the atoms (see Ref. [2J] for the pioneering 
use of this configuration). The two polarization modes 
are degenerate (simultaneously resonant in the cavity), 
so we can use the same cavity for both excitation and 
signal collection. 

We consider two different configurations for excitation. 
The first has a weak magnetic field parallel to the in- 
cident polarization. The light drives ir (Am=0) tran- 
sitions in the atoms in this basis (see Fig. [5^). An ex- 
cited atom can return to its original ground state m-level 



through spontaneous or stimulated emission into the par- 
allel (driven) mode; however, there is also some proba- 
bility to relax into a different m-level through a a + or 
<j~ (Am = ±1) spontaneous emission transition, which 
produces light that is circularly polarized with respect to 
the magnetic field axis and may be collected by the per- 
pendicular (undriven) cavity mode. Although the cavity 
coupling strength is modest, the amount of collected light 
increases measurably from that in free space [2|| . 

The second configuration has a weak magnetic field 
along the cavity axis, while keeping the same polariza- 
tion arrangement (see Fig. The linearly polarized 
drive can then be considered as the sum of two circu- 
lar components with different indices of refraction due 
to the Zeeman shift. This results in a Faraday rotation 
of the drive into the perpendicular mode. Though ordi- 
narily the Faraday rotation imparted by a single atom 
in a weak magnetic field is vanishingly small, the cavity 
finesse allows the light to compound multiple rotations 
before escaping, and the effect becomes appreciable. 

A polarizing beam splitter (PBS) sends the two modes 
to separate detectors, which can resolve individual pho- 
ton events. From the record of photon detections, we 
analyze the statistics of emission to look for signs of sin- 
gle atom transits. One advantage of collection into a 
cavity mode is that we can spatially filter most sources 
of background light. However, there are remaining con- 
tributions to the perpendicular mode from cavity mirror 
birefringence and background light that cannot be com- 
pletely eliminated. The optimal parameters for operation 
minimize the background influence on the atom detection 
confidence. 



III. THEORY 

The coupling of a single photon to a single atom in 
optical cavity QED takes place through electric dipole 
transitions at the Rabi frequency go = p, ■ E/h, where \x 
is the atomic dipole moment and E = y/hLo/2eoV x is the 
electric field amplitude of a single photon in the cavity 
mode volume yand polarization x. The coupling con- 
stant g varies in space with the cavity mode function. For 
a TEMoo Gaussian standing-wave mode of wavelength A 
in a cavity with axis along the z direction, the coupling 
takes the form: 

g(x 1 y 1 z) = g Q cos(kz)e-( x2+ y 2 y w K (1) 

where the amplitude go depends on the transition dipole 
moment and Clebsch-Gordan coefficients, k = 2n/X and 
wo is the mode waist. 

The real structure of the atomic hyperfine levels in 
our experiment is complex (Fig. [3^) with a wide range 
of values for go- We present here a three- level model 
(Fig. [3Jd) following Ref. [25|: two ground states (|1) and 
|3)) and one excited state (|2)). We call the coupling 
constants g and G for the two modes and the decay rates 
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FIG. 3: a) Diagram of 85 Rb levels relevant to the detection, 
b) Energy level diagram of the simplified model. 



7 and T for the two channels. The total decay rate of the 
population inversion is the sum, 



itot = 7 + r. 



(2) 



The rates are related by g/G = j/T, where the ratio de- 
pends on the Clebsch-Gordan coefficients for the specific 
transitions involved. The number of photons needed to 
saturate the driven atomic transition is: 



»o 



Itot 

33o 2 



(3) 



which includes a geometric factor from the Gaussian 
standing wave [26j . 

Together with the saturation photon number no, we 
take as a figure of merit the single-atom cooperativity 
C\ 1 a dimensionlcss parameter which characterizes the 
influence of a single atom in the driven mode with decay 
rate n: 



Ci 



An equivalent parameter for the undriven mode is: 



Cr = 



r* 2 



(4) 



(•5) 



We define a total cooperativity parameter C for N atoms 
in the two-mode system, which is slightly different from 
that for two-level atoms (i.e. C\N) [271 ] : 



C = 



dN 
l + 2Ci 



(6) 



We consider a cavity resonant with both atoms and 
input drive, and N atoms stationary and maximally- 
coupled to the mode. For a driving intensity of Y = 
(n) /no (representing the mean steady-state number of 



photons inside the empty cavity, normalized to the satu- 
ration number) and the three-level atom, two-mode cav- 
ity model discussed above, the steady-state photon num- 
ber with atoms in the driven (X\\) and perpendicular 
(Xx) modes is [251 ] : 



X n = Y 



1 



X\ 



Y 



(1 + 2C) 2 

2Ci 
l + 2Ci 



r(i-4C), 
c 

[l + 2Cf 



(7) 

Y2C 1 C, (8) 



where the lowest order approximations are valid for 2C -C 
1 and 2Ci -C 1. The signal X± increases linearly with 
the number of atoms and contains no additional terms, 
indicating that any amount of light is the result of atomic 
emission. This is in contrast to the parallel mode signal 
X|| which is the difference of the drive, determining the 
shot noise, and a small number proportional to 2C that 
carries the information of the atoms. 

To include the effects of Faraday rotation, we consider 
an atom with one excited state (m' = 0) and three ground 
states (m = ±1,0). (See Fig [2)3). A magnetic field B 
along the direction of light propagation lifts the ground 
state degeneracy through the Zeeman effect, resulting in 
a shift in the dispersion curves and a difference in the 
phase velocities of the two circular components of light, 
causing Faraday rotation [28| . We excite the system with 
vertically polarized light, which is an equal superposition 
of right- and left-circular ly polarized light. The light in- 
teracts with the atom and becomes elliptically polarized; 
the plane of polarization rotates through an angle. The 
light then passes again through the atom after reflection 
on the mirror and compounds the rotation. When the 
light exits the cavity, the total rotated angle is that of a 
single pass multiplied by the number of times the light 
passes through the atom, indicated here by the presence 
of the cooperativity C: 



2g L u B B/h-y t , 



1 + {2g L n B B/hj t , 



(9) 



where <?l is the Lande factor and \xb is the Bohr magne- 
ton. 

The addition of this effect changes Eq. [8] into: 



X± =Y (2C 1 C+\(j)\ : 



(10) 



Equation [TO] shows that the signal in the perpendicu- 
lar mode now consists of two parts: a background level 
of spontaneous emission (which may be considered as 
(cr + — er~) polarization if the drive is (a + + er~), in full 
analogy to the basis with tt drive), and Faraday- rotated 
drive which increases from zero to a maximum value of 
^±Faraday = Y C 2 / '4 when B = Wjtot / {^Ql^b) ■ For sim- 
plicity in what follows we treat the atom-cavity coupling 
in the basis of 7r drive, which must agree with the total 
coupling in the basis of a drive when the ground-state 
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populations are appropriately transformed (though the 
values of go and C\ for the individual transitions are gen- 
erally different). 

We analyze the two-time intensity autocorrelation 
function of the light from the perpendicular mode to 
understand the dynamics of photon emissions that con- 
tribute to the steady-state value of Xj_. We do not dis- 
tinguish between the two components in Eq. I10| treating 
both as arising from near-resonant scattering of incident 
light (resonance fluorescence) in which the cavity serves 
only as a mode in which to collect the light. 

For a stream of photons scattered from a single two- 
level atom undergoing stationary, resonant excitation, 
the intensity autocorrelation function takes the familiar 
form for resonance fluorescence antibunching: 

<7^ 2) (t) = 1 - e^/ 4 ^ (cosher + ^ sinhSr), (11) 



9 {2) {r) = l + 



f{r)9 ( 2\r) 



where 8 = (7/4) vT - -~ 8Y is a function of the drive inten- 
sity Y 29]. The antibunching referred to here means an 
increasing slope for increasing times: dg^(r) / 'dr| T >o > 
0. The condition g^ 2 \0) < 1 indicates that the photons 
obey sub-Poisson statistics, and is not a necessary con- 
dition for the photons to exhibit antibunching. 

The modification of Eq. [11] for the case of an atomic 
beam passing through a localized excitation region ap- 
pears in Ref. [HJ, which adds a number of effects: L, the 
multiplication of Eq. [11] by a window function to include 
the shape of the mode function encountered by each tran- 
siting atom; ii., the scaling of Eq. [TT]by the inverse of the 
steady-state mean number of interacting atoms, N; iii., 
the addition of a term \gy (r)| 2 to include the beating 
of fields emitted by different atoms; iv., the addition of a 
constant offset of one arising from uncorrelated emissions 
from different atoms; and v., the reduction of the correla- 
tion factor due to the presence of completely uncorrelated 
background light. Further generalizations appear in Ref. 
[22I . I23I ] , where the authors also consider the product of 
the fields that come from atomic emission with a tempo- 
rally uncorrelated but spatially mode-matched field from 
background light. 

Our system is more complicated than either of these 
cases because the birefringence background, as a compo- 
nent of the driving laser, is temporally coherent and com- 
pletely mode-matched with the atomic emission. This 
gives a non-zero expectation value for the product of 
atomic and background fields, and adds a number of 
complicated terms to the correlation function which be- 
come increasingly important with increasing drive inten- 
sity. Such terms should decay to zero on the time-scale 
of the atomic lifetime, however, and are unimportant for 
weak driving, so we omit them from the present analy- 
sis and concern ourselves with understanding the atomic 
motion and number fluctuations, which manifest them- 
selves mainly in the long tails of the correlation function. 

We are able to describe the total correlation function 



(l + R b /R s ) 2 N 



-F{T;R b /R s ,N) 
(12) 

where R b and R s are the average count rates for back- 
ground and signal, respectively, and 



/(r) = [cos(fir)e-^ T + \]e^ T ' T ^ 



(13) 



is an empirically-determined window function describing 
motion of the atom through the cavity mode with Gaus- 
sian temporal width T, where a small tilt of the atomic 
beam with respect to the cavity axis normal results in 
motion through the standing-wave lobes at frequency O. 
These oscillations decay at a rate f3 much faster than 
the transit time, due to the spread in velocities from the 
atomic source. The function F corresponds to any ad- 
ditional terms related to the beating of signal from dif- 
ferent atoms or signal with background, which at low 
intensities is dominated by partially mode-matched scat- 
tered light from outside the cavity. Since this term repre- 
sents classical noise sources, it can only contribute light 
which is bunched around r = 0. Any observation of an- 
tibunching in our signal must come from the single-atom 
(2) 

term g\ ' (r) . An explicit observation of F would require 
N ^S> 1 since its width is comparable to that of the anti- 
bunching signal [181 ]. 



IV. APPARATUS 

The apparatus consists of two main components: the 
source of atoms and the cavity. A titanium sapphire laser 
(Ti:Sapph) provides most of the light needed for the ex- 
periment at 780 nm. The laser linewidth and long-term 
lock are controlled using a Pound-Drever-Hall (PDH) 
technique on saturation spectroscopy of 85 Rb. The ex- 
citation beam of the cavity passes through an electro- 
optical modulator to imprint frequency sidebands sepa- 
rated by 150 MHz from the carrier. One sideband is near 
resonance with the atomic transition and the cavity, be- 
coming the system drive, while the carrier and the other 
sideband are far detuned and reflect from the front cav- 
ity mirror. Changing the sideband frequency allows us to 
probe the frequency response of the resonant cavity-atom 
system without misalignment of the input laser coupling. 

Polarization elements (better than 5 x 10 -5 ) and mode- 
matching optics prepare the driving laser before it enters 
the cavity. A lens at the exit collimates the beam and a 
half-wave plate (HWP) aligns the polarization to a calcite 
PBS which separates the parallel (driven) and perpendic- 
ular (undriven) modes. The perpendicular mode passes 
through a second beam splitter which can divide the 
light between two avalanche photodiodes (APD) (Perkin 
Elmer SPCM-AQR-12 and -13) for single or coincidence 
measurements. A series of filters, telescopes and aper- 
tures before the detectors remove background light. The 
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APD electronic output pulses then go to a time-stamp 
unit in a computer or to a correlator (30j . The system 
has a 23% measured photon detection efficiency for the 
APDs and optical paths. 

The Fabry-Perot cavity has a mirror separation of 2.2 
mm and a 1/e field TEMoo mode waist of 56 /im. The 
input mirror transmission (15 ppm) is lower than the 
output mirror transmission (300 ppm) by a factor of 20 to 
ensure that most of the signal escapes from the cavity on 
the detector side. The decay rate for the cavity is k/2it = 
3.2 x 10 6 s- 1 with finesse of 11,000. The cavity length is 
kept resonant with the F — 3 — > F 1 — 4 transition of the 
D 2 line of 85 Rb at 780 nm. The birefringence splitting 
of the two polarization modes is less than 500 kHz. The 
cavity length is stabilized by the PDH method with light 
derived from an auxiliary laser at 820 nm. A grating and 
interference filters separate the 780 nm signal light and 
820 nm locking light on the cavity output. 

A rubidium dispenser delivers Rb vapor to a magneto- 
optical trap (MOT) in a stainless steel chamber above the 
vacuum chamber that houses the cavity. Both are kept 
under vacuum with pressures below 10 -8 Torr. We use 
a six-beam configuration with 1/e power diameter of 20 
mm and 10 mW per beam. A second laser repumps the 
atoms that fall out of the cycling transition in the trap. 
Combinations of acousto-optical modulators (AOM) per- 
mit independent settings of the frequency and amplitude 
of all laser beams, controlled by computer. A pair of 
coils generates a magnetic field gradient of 10 G/cm, and 
three sets of independent coils zero the magnetic field at 
the trapping region. Three extra coils around the cavity 
region allow for configuring a constant field perpendicu- 
lar or parallel to the axis of the cavity without disturbing 
the location of the MOT. 

The retro-optics of the vertical arm of the MOT (in- 
side the vacuum chamber) have 1.5 mm-diameter holes 
on axis. This creates an imbalance in the light pressure 
of the vertical cooling beam such that cold atoms are 
pushed down into the cavity, producing a low velocity 
atomic beam (LVIS) [3l|. The atomic beam intersects 
the cavity mode about 8 cm below the trapping region. 
The mean speed of the atoms passing through the cavity 
mode is approximately 15 m/s. Ref. [3l[ shows that the 
distribution of longitudinal speeds is considerably wider 
than that of a sample at the Doppler temperature, while 
the spread in transverse velocities results primarily from 
the geometric collimation. 

The single-atom dipole coupling frequency for the 
driven cavity mode depends on the Clebsch-Gordan coef- 
ficient for the particular states involved. Table 1 gives the 
expected range of values for go, C\, and no, for transitions 
between the F = 3 ground state and the F' = 4 excited 
state. The strongest coupling for 7r-polarized light occurs 
for to = and the weakest for to = 3. The atoms enter 
the cavity in a distribution of ground states that is not 
centered around m = 0. The total atomic decay rate is 
7 tot /27r = 6 x 10 6 s _1 . We write all intensities normalized 
to the m = 0, 7r-polarized saturation intensity, hq = 5.3 



TABLE I: Clebsch-Gordan (CG) coefficients, single-atom cou- 
pling constants go, single-atom cooperativities C\, and sat- 
uration photon numbers no for different transitions in the 
D 2 line of 85 Rb from F = 3 to F' = 4, for different drive 
polarizations(7r or tx + ). 





F = 


3 


F' 


= 4 


CG 


go/2-K [MHz] 


Ci 


n 




m — 





m' 


= 




1.5 


0.12 


5.3 




m — 


3 


m! 


= 3 




0.99 


0.053 


12 




m — 





m! 


= 1 


v/5/28 


1.2 


0.075 


8.5 


a+ 


m — 


3 


m! 


= 4 




2.0 
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-3 
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Vl/56 


0.38 
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photons. 

We find the peak count rate for the perpendicular mode 
with Faraday rotation by scanning the laser and cavity 
together across the atomic resonance. We take all mea- 
surements at the location of the peak (|S|=3.3 Gauss), 
which is shifted by about 3 MHz from the F = 3, to = 
— > F' = 4, to' = transition. 

For evaluating photon statistics, we detect the perpen- 
dicular mode light on a single APD for 300 s at different 
drive intensities. The detection records are stored on 
computer for post-processing and analysis, with a reso- 
lution of 4 ps. For measuring the intensity autocorrela- 
tion, we split the light between two APDs to eliminate 
distortions from detector dead time and after-pulsing. 
These are also recorded for 300 s at each intensity, at 
the same atomic beam density used for the single APD 
measurements. The autocorrelation is formed by mak- 
ing a histogram of the time between coincidences on the 
two detectors, out to a delay of ±10 us with a bin res- 
olution of 10 ns. We take as normalization the mean 
number of counts per bin, based on the long-time aver- 
age count rates. Measurements for calibrating Faraday 
rotation were taken separately, at a similar atomic beam 
density. 



V. PRELIMINARY DATA ANALYSIS 

A. Faraday rotation 

The Stokes parameter formalism [32j allows a measure- 
ment of the rotation angle <j> if the amplitudes of the elec- 
tric field in the two polarizations e±_ and eu are known: 

w = F^HTP * r • (14) 
|eiir-|e±r £ \\ 

where the approximation is made for |eii| |ej_|< We 
measure the count rates in the perpendicular mode as a 
function of applied magnetic field along the cavity axis 
to determine the location of maximum Faraday rotation. 
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The peak rotation is different depending on the sign of 
the magnetic field, as the atoms experience different op- 
tical pumping by the pushing beam from the MOT, and 
we do not otherwise optically pump before the cavity. 
We obtain a maximum of \cj>\ — 0.035 ± 0.007 rad from 
the measured increase in the count rate with a field mag- 
nitude of 3.3 Gauss. To make explicit the effects of the 
cavity, we estimate the effective number of maximally- 
coupled atoms by using Eq. [7] and measuring the driven 
mode transmission. This was greater than 0.95, giving a 
bound of C <0.01 or N <0.1 effective maximally-coupled 
atoms (for ir transitions from to — 0). Calculating the 
free-space absorption length for 0.1 atoms and our cav- 
ity waist, we would expect a maximum rotation of the 
order of 3 /xrad [28(. The enhancement by a factor of 
approximately 10,000 is a result of the cavity finesse. 



B. Photons per atom 

The formula of Mandel relates the photon number 
distribution P(n) with the atom number distribution 

Patom{m): 




^0.4 



1.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

<n> 




P{n) = P(n\m)P at 

Oil! 

(m), 



(15) 



FIG. 4: Examples of data (circles) and least-squares fits (con- 
tinuous line) to Eq. (|16|l for Faraday rotation (top) and for no 
Faraday rotation (bottom) as a function of average number 
of photons in a time bin. 



where P(n\m) = (am)™ exp (— am)/n\ is the conditional 
probability of detecting n photons when there are m 
atoms in the cavity volume, each contributing a mean 
of a photons to the signal with a Poisson distribution of 
number (33|. This shows that the super-Poisson fluctu- 
ations of the light arise from the combined Poisson fluc- 
tuations in atom number and photons scattered, with a 
an extractable parameter. 

We follow the method of Ref. [II] to relate the mean 
and mean squared photon number ((n), (n 2 )) to a: 



1 = {n)9aa + OL, 



(16) 



where the atom-atom correlation function is g aa = 
((to 2 ) — (to)) /(to) 2 , which for Poisson atomic fluctua- 
tions reduces to g aa = L 

Figure Q] shows a detail of the plot of the experimen- 
tal values for the left side of equation p^|) with the two 
different detection methods: Faraday rotation on the top 
and no Faraday rotation on the bottom. We calculate 
the statistics based on time series of 300 s that we bin 
in time. Longer bin sizes have a larger average of counts 
(horizontal axis). The intercept gives the value of a. A 
linear fit for bins from 50 us to 100 us ((n) = 1.81 to 
3.63, off the scale of the plot in FigHJ) in the Faraday 
case gives a vertical axis intercept of a = 0.196 ± 0.003 
photons detected per atom, and a slope of 1.033 ± 0.001. 
The intercept in the no Faraday case (for a fit from (n) = 
0.50 to 1.00) is a = 0.036 ± 0.002 photons detected per 
atom and a slope of 1.026± 0.003. The drive in both 



cases is 0.4 no. The value of the slopes indicates that 
the assumption that the atoms follow Poisson statistics 
is well justified. 

Figure [5] shows the extracted values of a as a func- 
tion of cavity photon number (driving intensity) span- 
ning more than two orders of magnitude. There are 
data points (rhombs) that come from Faraday rotation 
and one with spontaneous emission alone (circles). The 
latter is less efficient by almost a factor of five, a in- 
creases linearly with drive until atomic saturation inten- 
sity (n ~ no), where also optical pumping effects start to 
enter. Each point on Fig. [5] comes from a least-squares 
fit similar to those in Fig. 2] The count rates in the per- 
pendicular mode and the values of a give an atomic flux 
of approximately 160,000 atoms s _1 . 



VI. ANALYSIS AND OPTIMIZATION OF 
DETECTION 

Three identifiable sources of background counts con- 
tribute to the measurement of the perpendicular mode: 
intrinsic detector dark counts (approximately 300 s _1 
each), scattered light from the MOT beams (approxi- 
mately 2000 s _1 ), and the birefringence of the cavity 
(less than 5 x 10~ 5 ). The light from the MOT dominates 
the background count rate at low driving intensity; how- 
ever, at higher intensity, light from cavity birefringence 
is the main contributor to the degradation of the signal- 
to-background ratio. A single photon escaping in the 
perpendicular mode most likely comes from an atom in 
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FIG. 5: Counts detected from a single atom (rhombs for Fara- 
day, circle for non Faraday) as a function of intensity of the 
exciting laser normalized by the saturation photon number. 
The arrows indicate the points extracted from the plots on 
Fig. H] The error bars are smaller than the size of the data 
points. 
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FIG. 6: Signal-to-background ratio as a function of normal- 
ized driving laser intensity for single detection (open rhombs) , 
double coincidence detection (open squares), and triple co- 
incidence detection (open triangles) in a 1 /is window with 
Faraday rotation of the drive. 



the cavity; however, there is at best still a 4% probability 
that it comes from background counts. We suppress the 
probability of a false atom detection by requiring photon 
coincidences in a time window smaller than the transit 
time of an atom across the cavity mode 16]. It is pos- 
sible to implement this in real time using an electronic 
coincidence counter. The first perpendicularly polarized 
photon detection ( "start" ) opens a gate of variable width 
which allows counting pulses from the "stop" APD. 

Figure [6] shows the signal-to-background ratio calcu- 
lated from time series (a posteriori in this case) for 
the rate of single photon detection (open rhombs), two- 
photon coincidence (open squares), and three-or-more- 
photon coincidence events (open triangles) in a 1 [is win- 
dow as a function of driving intensity. The maximum 
occurs at the same driving intensity for all three. Two 
photons within a 1 /is window improves the signal-to- 
background ratio by more than an order of magnitude 
compared to detection of single photons. The three- 
photon coincidence gives a better ratio but significantly 
decreases the rate of detection as exemplified by the error 
bars associated with each point. The three sets achieve 
their maxima just before the atomic saturation inten- 
sity, due to background counts from cavity birefringence 
that continue to increase linearly with drive. The contri- 
bution to the coincidences from detector afterpulsing is 
small (less than 1%) and does not affect the results. 

To evaluate the probability of making a false detec- 
tion for various coincidence gate times, we extract the 
fidelity of detection F versus gate length for two-photon 
coincidence based on the time record. We calculate the 
waiting time distribution of photon arrivals on a single 
APD by making a histogram of the time separation be- 
tween consecutive detections. The integral of the waiting 
time distribution gives the number of two-photon coin- 



cidences c in a given time window. We calculate c with 
and without atomic flux to obtain: 

F 1 (-without 1 1 ^with- (^-^) 

The fidelity (one minus the probability of error) 
reaches more than 99% at 0.1 /is, and 99.7% at 1 ^s with 
Faraday rotation, but only 96.7% at 1 fis with detection 
of spontaneous emission alone. The fidelity decreases at 
longer times. The fidelity is optimal for times between 1 
and 5 /xs, determined by the distribution of atom transit 
times. 



VII. INTENSITY AUTOCORRELATION 

We ensure that our coincidence detection scheme is 
sensitive to single atoms by measuring the intensity auto- 
correlation function (g^( r )) of the perpendicular mode 
under very weak driving intensity at the same atomic 
beam density used during detection measurements. If we 
consider the light in the perpendicular mode as coming 
from resonance fluorescence, the antibunching can only 
arise from photon pairs from the same atom. The obser- 
vation of antibunching in our signal is strong indication 
of our ability to detect photon coincidences from a single 
atom in less than 1 /zs. 

Figure [7] shows two examples of an autocorrelation 
function around r = with and without Faraday rotation 
for the same driving intensity (0.4 uq). The antibunching 
is visible and lasts for a time of the order of the excited 
state lifetime (26 ns). The area under the curve is larger 
with Faraday rotation, indicating a substantial increase 
in photon flux from individual atoms. The antibunching 
is less pronounced without the rotation because of the 
lower signal-to-noise ratio. 
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FIG. 7: (Color online) Intensity autocorrelation function 
(g (t)) of the perpendicular mode with (triangles) and with- 
out (diamonds) Faraday rotation, at a driving intensity of 0.4 
no. 
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FIG. 8: (Color online) Intensity autocorrelation function 
(g (t)) of the perpendicular mode at 0.24 no. The con- 
tinuous line is the theoretical fit to Eq. I12I 



Since we record the photon detections as continuous 
time series, we can also look at the long-term behavior of 
the autocorrelation in the context of resonance fluores- 
cence as presented in Ref. [U [22], [23} . We fit the data to 
Eq. [12] with free parameters N, T, 0, and f2, excluding 
the term F(t) and the region |r| <50 ns where it is im- 
portant. We find excellent agreement in this outer region 
(reduced x 2 =l-03 for points out to ±5 fis). Fig. [5]shows 
the fit to the data taken at a driving intensity of 0.24 n sa t- 
The damped oscillation at f2=1.5 MHz corresponds to 
motion through the standing wave with a period of 0.67 
fis (a velocity of 0.58 m/s along the cavity axis) and a 
damping time of l//3=0.29 /is. Monte Carlo simulations 
of the transit show good agreement with this damping 
time based on the geometric collimation of transverse ve- 
locities from the source [3l|. The Gaussian background 
gives a mean atom number of iV=0.88, and a 1/e waist of 
T=2.7 /is (a mean velocity of 14.7 m/s across the Gaus- 
sian mode) . The ratio of the velocities gives a beam tilt 
of 2.3° with respect to the cavity axis normal. Although 
we take care to align the cavity under the exit hole of the 
atomic source, the cavity mode is not necessarily cen- 
tered on the two mirrors, and the transverse beam width 
allows for small inclinations of the beam. 

The full structure of the correlation function depends 
not only on the above parameters (properties of the 
atomic beam), but also on the intensity of the drive, 
which can change the atomic response or increase the 
amount of background light from cavity birefringence. 
Fig. [5] shows the evolution of the autocorrelation func- 
tion as the driving intensity increases by two orders of 
magnitude. We observe a transition from antibunching 
to bunching with higher drive. The behavior comes from 
the term F(t) in Eq. [T2] which includes beating of the 
birefringence light with the atomic emissions. This is 
in contrast to the transition seen in Ref. [ll[, where 
the atomic density is increased such that the constant 
\9a ( t )| 2 t erm (part of our F(t)) becomes visible. On 
separate experiments at low intensity we have also fol- 




FIG. 9: Intensity autocorrelation function (gr 2 - 1 (r)) of the per- 
pendicular mode for different values of photon number (inten- 
sity) in the driven mode. 



lowed the disappearance of the antibunching as we in- 
crease the number of atoms. 

Although the actual density of atoms is constant, 
the effective excitation volume increases with power, as 
atoms near the nodes and wings of the mode start to 
interact [2^|. This process reaches a maximum above 
saturation intensity, when the only unsaturated atoms 
are too weakly coupled to emit into the mode. Fig. 
[TDk . shows the extracted values of TV versus driving 
intensity for the autocorrelations in Fig. [9] The line 
shows the expected shape of the saturation, assuming 
that the probability of collecting a photon from an atom 
is proportional to the product of a constant g 2 {x,y,z) 
for collection and a saturating strength of excitation 
(i(n/n )g 2 (x, y,z)]/[l + (n / 'n )g 2 (x , y , z)]), which is con- 
sistent with the expression Y2C\C in Eq. [5] modified for 
strong driving. We obtain the theoretical line by inte- 
grating this product over a volume much larger than the 
cavity mode, and using a least-squares fit to scale the 
vertical axis and no. 

Using the extracted parameters of mean atom num- 
ber and mean transit time from the two-APD correla- 
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FIG. 10: a) Measured mean number of interacting atoms, 
and theory curve showing the expected saturation, b) Mea- 
sured (filled rhombs) count rates in the perpendicular mode 
together with those predicted from the extracted mean num- 
ber of atoms (open squares), photons per atom, and transit 
time. 



tion measurements, together with the extracted values 
of a from the one-APD measurements, we can predict 
the expected macroscopic count rates as R s = Na/2T. 
Fig. [10b . shows that we obtain excellent agreement for 
low intensities, while the highest intensity drive (largest 
birefringence background) disagrees significantly. This 
comes from the large contribution from the beating be- 
tween background and signal, which generates the large 
central peak in the measured autocorrelation. 

We emphasize that although the mean number of inter- 
acting atoms is approximately one, this is not the same as 
the effective number of maximally-coupled atoms, since 
most are weakly coupled. To explore this numerically 
we use the work of Carmichael and Sanders (34[ • We dis- 
tribute atoms randomly and uniformly across the cavity 
mode function such that the mean density is 0.88 atoms 
within two mode waists (similar to the measured value 
for low driving intensity). Allowing the atom number to 
fluctuate with a Poisson distribution, we sum the indi- 
vidual coupling strengths for each configuration of atoms 
and weight by the atomic number distribution to obtain 
the mean steady-state effective number of atoms, which 
in this case is 0.04. Using the mean value of Ci=0.09 for 



7T transitions from the ground state, this gives a driven 
mode transmission of approximately 0.99 from Eq. 
consistent with the experimental conditions for single 
atom detection. A rapid decrease in probability above 
one effective atom indicates that we operate in a density 
regime of single-atom coupling, with a low probability of 
having two atoms simultaneously coupled by more than 
half of the maximum value. This numerical result to- 
gether with the antibunching in the perpendicular mode 
counts confirms that the measured signals indeed arise 
primarily from single-atom emission bursts. 



VIII. CONCLUSION 

A two-mode cavity QED system is an excellent single- 
atom detector. The detection of a single photon out of 
the undriven mode signals the arrival of an atom in the 
cavity, and the detection of a second photon within 1 /is 
confirms the presence of an atom with more than 99% 
confidence. The method can operate over a wide range 
of photon fluxes and works for atoms traveling as fast as 
15 m/s across the mode. Because the cavity is frequency- 
selective and has a short lifetime (approximately equal to 
the atomic excited state lifetime) , the method guarantees 
that the atom is in the F — 3 hyperfine ground state 
immediately after the detection of the second photon. 
In principle information on the m-state could be gained 
through initial optical pumping and increased Zeeman 
separation of levels to distinguish the transitions in fre- 
quency. However, this would also require clearly choosing 
a detection method between Faraday rotation (no change 
in m) and spontaneous emission collection (Am = ±1). 
In the regime where we have demonstrated detection, the 
combination of the two processes leaves the atom in a 
mixed state of m-sublevels after the photon detections. 
This technique is a classical determination of atom pres- 
ence rather than a projection onto a pure quantum state. 

One significant advantage of this scheme is that it is 
sensitive enough to work without strong coupling of the 
cavity and atoms, meaning that the detection volume 
can be kept large. As presented in Eq. [101 the signal 
rate for both spontaneous emission collection and Fara- 
day rotation (for small rotations) depends on the square 
of the single-atom cooperativity, which is proportional 
to the cavity finesse. For a given mode volume, increas- 
ing the finesse by some factor (through higher reflection 
coefficients on the mirrors) will in general increase the 
detection efficiency. However, the cavity linewidth de- 
creases by the same factor, which may reduce the rate of 
escape of light from the cavity. The signal-to-noise ratio 
of Faraday rotation to birefringence is independent of the 
finesse, since both arise from a rotation of the polariza- 
tion on each round trip of the light. 

We have demonstrated how autocorrelation measure- 
ments of the emitted light show the underlying dynamics 
of the atom as it traverses the Gaussian mode of the 
cavity and encounters the standing wave. This method 
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of atomic detection, using a conditional measurement, 
is well-suited for some quantum control protocols. The 
implementation should facilitate the manipulation of the 
electromagnetic field produced in the atom-mode inter- 
action of cavity QED systems. 
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